The International Journal Of Engineering And Science (IJES) 
II Volume 113 II Issue II 4 II Pages II 65-70 II 2014 II 
ISSN (e): 2319 - 1813 ISSN (p): 2319 - 1805 

Investigation of behaviour of thermodynamic properties of 
cassava based ethanol-water system during azeotropic distillation 

process. 

lj A. E. Adeleke, 2 P. O. Aiyedun, 3 O.U. Dairo 

' Mechanical Engineering Department, Osun State University, Osogbo, Nigeria. 
2 ' Mechanical Engineering Department, Federal University of Agriculture, Abeokuta, Nigeria 
' Agricultural Engineering Department, Federal University of Agriculture, Abeokuta, Nigeria. 

ABSTRACT 

This paper investigates the behaviour of properties of cassava based ethanol-water system during azeotropic 
distillation process. The product of fermentation is assumed to be binary mixture. The objective is to arrive at 
the optimum process conditions that will give the required percentage of purification using a new developed 
software. Simulations are performed with cassava as the biomass material and the parameters obtained from 
the beer produced from fermentation of cassava mash as the input for the software.The considered parameters 
were Reflux Ratio (R), Relative Volatility (RV), Minimum Reflux Ratio (R min ), Mole Fraction of the most volatile 
component in the distillate (x D ), Mole Fraction of the most volatile component in the bottom (x B ), Mole Fraction 
of the most volatile component in the feed (x F ), thermodynamics properties of the feed (q). The results obtained 
using the developed Software is compared with the predicted result from Central Composite Design model at 
various input levels and the new software is validated with published experimental results. The process 
simulation results show that the reflux ratio, relative volatility increases with increase in temperature at 
constant pressure and have significant effect on purification of ethanol compared to other parameters. 
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I. INTRODUCTION 

The search of environment friendly solutions for the management, along with increasing costs of fuel 
and recent regulations on disposal, emissions leads toward the recovery of energy and demands research 
activities related to renewable energy such as biofuel which comprise of bioethanol and biodiesel [1]. Bio- 
ethanol is a distilled liquid produced by fermenting sugars from sugar plants and cereal crops such as sugarcane, 
corn, beet, cassava, wheat and sorghum. Several researchers have reported biofuel production from various 
agricultural materials including biogas from mixture of cassava peels and livestock waste [2], fuel from 
indigenous biomass wastes [3], as well as ethanol from non-edible plant parts [4]. 

Ethanol is a common ingredient in many solvents, perfume, disinfectants and alcoholic beverages. To 
produce bioethanol, it is necessary to separate ethanol from ethanol-water mixture after fermentation of alcohol. 
An azeotropic distillation is a widely practiced process to split one needed component from its heterogeneous 
system where the third component known as entrainer is incorporated. Azeotropic distillation columns are 
commonly used in industry to separate mixtures of close relative volatility and for breaking azeotropes such as 
bioethanol -water system. Process simulation is used for the design, development, analysis and optimization of 
the technical processes. It is a model-based representation of chemical, physical, biological, energy producing 
and transforming, in addition to other technical processes and unit operations using software. Process simulation 
always uses models which introduce approximations and assumptions, but allows the description of a property 
over a wide range of parameters which may not be covered by real data. Models also allow for interpolation or 
extrapolation, searching for conditions outside the range of known properties [5,6]. Many researchers have 
studied the production of bio ethanol from different feedstock but not much have taking coginizance of the 
behaviour of thermodynamic properties of the feed during distillation process. This paper investigates the 
behaviour of thermodynamic properties of cassava based ethanol-water system during distillation process. 
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II. MATERIALS AND METHODS 

A new developed software [7] hinged on McCabe-Thiele method and the Fenske Equations [8] were 
used for the simulation of thermodynamics properties of azeotropic distillation of cassava based bioethanol 
method using Vapour- Liquid- Equilibrium (VLE) data of ethanol-water system obtained from literature [9]. 
The screen shot of the software used is as shown in Figure 3. Table 2 shows the output values obtained 
from.The Design Expert software was used to obtain the output values at different levels in the developed 
software. The central composite design is used to compare the values obtained from the developed software and 
the predicted values. 




Figure 1: Screen shot of Software 



The model was used by supplying the input and the level of simulation was provided which gives an 
output that translated to the plotting of the McCabe-Thiele graph.The report is generated and displayed. 
Parameters such as RR, Rc 0 eft> x d> x b> x b Viscosity, q were varied in the programme to obtain the Number of 
trays. The values obtained from the software were verified by comparing to values obtained from mathematical 
computation. To obtain the optimum design parameters for use in mechanical design, a central composite 
experimental design was used. 

McCabe-Thiele Theoretical Analysis 

The McCabe-Thiele method uses the equilibrium curve diagram to determine the number of theoretical stages 
(trays) required to achieve a desired degree of separation. It is a very useful tool for the understanding of distillation 
operation [8] .The VLE data is available at the operating pressure of the column. Information required were the feed 
condition (temperature, composition), distillate and bottom compositions; and the reflux ratio (R) is defined as the ratio of 
reflux liquid over the distillate product. The feed has a concentration of x F (mole fraction) of the more volatile component, 
and a distillate having a concentration of X D of the more volatile component and a bottom having a concentration of x B was 
desired.The method involves the plotting on the equilibrium diagram, 3 straight lines: the rectifying section operating line 
(ROL), the feed line (also known as the q-line) and the stripping section operating line (SOL). Each of these lines passes 
through the points representing the mole fractions of the more volatile component in the distillate, bottoms and feed (x D , x B 
and x F ) respectively. These lines represent the relationship between the concentrations in the vapour phase (y) and the liquid 
phase (x). The number of theoretical stages required for a given separation was then the number of triangles that can be 
drawn between these operating lines and the equilibrium curve. The last triangle on the diagram represents the reboiler. 



In the simplest case, the McCabe-Thiele Method to determine the number of theoretical stages followed the steps below: 

• Analysis of the Rectifying section, and determining the ROL using x D and R 

• Analysis of the Feed section, and determining the feed condition (q) 

• Determination of the feed line (q-line) using x F and q 

• Locate the intersection point between ROL and q-line 

• Analysis of the Stripping Section, and determining the SOL using, ROL, q line and x B 

The plates are numbered serially from the top and the plate under consideration is the n-th plate from the top. The 
plate immediately above this plate is plate (n-1) and the plate immediately below this plate is plate (n+1). 
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The Operating Line Equation for the rectifying section under constant molal flow is given as: 
R 1 

y " +l ~~{R^l) Xn+ ~{R^l) XD 1 

Usually the subscripts "n+1" and "n" are not known, so that equation 1 can be put together in the form equation 
2 

R 1 

y- x h x n 2 

(R + l) (R + l) 

The feed operating line equation is given as: 
q 1 

v = x-\ X F 3 

(1 + q) (l-q) 

The stripping operating line (SOL) equation is given as: 
L B 

y = x x R 4 

(L'-B) (L'-B) 

Total Reflux 

The Fenske Equation was also used to determine the minimum number of trays required for a given separation 
using eqn. 5. 



In 



X D (l X B ) 



ln(a AB ) 

Minimum Reflux 

The minimum reflux rate was determined mathematically from the endpoints of the rectifying line at 
minimum reflux — the overhead product composition point (x D ,x D ) and the point of intersection of the feed line 
and equilibrium curve (x', y'). 

K min ~ — 7 6 

y -x 

III. RESULTS AND DISCUSSION 

The results of the effect of Reflux Ratio and relative volatility on number of trays obtained using the 
developed software at constant temperature is shown in Figures 2 and 3 respectively. It was observed that 
increase in reflux ratio and relative volatility led to decrease in the value of number of trays. This is not 
unexpected as high refluxration during distillation process reqired less number of trays. [10] in a research on 
ethanol productivity from cassava reported similar observation in his study. The variation in relative volatility 
showed corresponding number of trays of distillation column that would give the highest purity of ethanol. The 
obtained values from the developed software and the predicted values from central composite design are also 
shown in Table 1. A correlation coefficient of 0.97 was obtained for the two results indicating that the 
developed software can adequately provide the expected output within the input range considered. A model was 
obtained relating the various experimental variables to the outputs; Number of Trays using Fenske and McCabe, 
as well as the Column height and Efficiency of distillation column. The model was then used to establish the 
optimum variable parameters for the distillation process. Model Optimization result obtained from the central 
composite design is shown in Table 2. 
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Table 1 . Comparison of results obtained from the Developed Software (DS) and predicted result from Central 

Composite Design (CCD) model at various input levels 
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Table 2. Optimization result obtained from the central composite design 
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Figure 2: Effect of Reflux Ratio on distillation column height 
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Figure 3: Effect of Relative Volatility on number of plates 

A partubation plot also revealed that the optimum values for parameters are close to the centre points. 
It also showed that R coe ff and RV are the two most important variables affecting our desire (Optimization), while 
Reflux ratio had the least influence (Figure 4). 
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Figure 4: Interactive effects Desirability and deviation from reference point on perturbation 

IV. CONCLUSION 

The behaviour of thermodynamic properties of ethanol-water system during the distillation process was 
studied and the parameters that will give the best grade of ethanol with minimal impurities were optimized for 
effective design of distillation column. The process simulation results show that the reflux ratio, relative 
volatility increases with increase in temperature at constant pressure and have significant effect on purification 
of ethanol compared to other parameters. 
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